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■ Abstract 

<D ■ 

■ We present a measurement of the strong coupling as using the three-jet rate 
measured with the Durham algorithm in e + e _ -annihilation using data of the JADE 
experiment at centre-of-mass energies between 14 and 44 GeV. Recent theoretical 
improvements provide predictions of the three-jet rate in e + e~-annihilation at next- 
to-next-to-leading order. In this paper a measurement of the three-jet rate is used to 
determine the strong coupling as from a comparison to next-to-next-to-leading order 
predictions matched with next-to-leading logarithmic approximations and yields a 

cj ■ value for the strong coupling 

a s (M z o) = 0.1199 ± 0.0010(stat.) ± 0.0021(exp.) ± 0.0054(had.) ± 0.0007(theo.) , 

co : 

■ consistent with the world average. 
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1 Introduction 



The annihilation of an electron-positron pair into a pair of quarks provides an ideal labora- 
tory to test the theory of the strong interaction, Quantum Chromodynamics (QCD) [2H3]. 
The free parameter of QCD, the strong coupling as, can be determined with events with 
more than two jets in the final state. To first order perturbation theory the radiation of 
a gluon from a quark is proportional to as- To determine this fundamental constant the 
observed three jet rate is compared to a perturbative expansion which predicts the three 
jet rate as a function of a single parameter as [SHI]- 

Recently theoretical progress has been made leading to a significant improvement 
in the prediction of the three jet rate (SHU] as a function of as- Previously e + e~-event 
shape distributions and the three-jet rate were only known to next-to-leading-order (NLO) 
accuracy, now QCD calculations to next-to-next-to-leading-order (NNLO) are available. 
These predictions were used to determine as based on a single value of the resolution 
parameter of the three-jet rate using the Durham algorithm [T2], with data taken at a 
centre-of-mass energy at 91 GeV with the ALEPH experiment [T3] . 

In this analysis we use data taken with the JADE experiment located at the PETRA 
collider at DESY between the years 1979 and 1986. We measure the strong coupling 
as at six centre-of-mass-energies in the range between 14 and 44 GeV. Besides using a 
different energy range we perform a fit in a range of the three-jet resolution parameter. We 
present a matching scheme to combine NNLO predictions together with next-to-leading 
logarithm approximations (NLLA). The matched predictions are used to determine the 
strong coupling as- The analysis follows closely the determination of as measuring the 
four-jet rate using data collected with the JADE experiment [T4] . 



2 Observable 

To determine a multijet-rate a jet finding algorithm has to be applied to particles observed 
in the final state. For this analysis we use the Durham algorithm [12], which selects jets 
according to the jet resolution parameter y cut . The Durham algorithm defines initially 
each particle as a proto-jet and a resolution variable yij is calculated for each pair of 
proto-jets i and j: 

2min(£2,£ 7 2 ) 

va = — ^2 ( x - cos %)> (!) 

Ads 

where Ei and Ej are the energies of jets % and j, cos 9ij is the cosine of the angle between 
them and E vis is the sum of the energies of the detected particles in the event (or the 
partons in a theoretical calculation). If the smallest value of ^ is less than a predefined 
value y cut , the pair is replaced by a new proto-jet with four-momentum p% = p^ + Pj, and 
the clustering starts again. Clustering ends when the smallest value of is larger than 
?/ cu t, and the remaining proto-jets are counted as final selected jets. 

Perturbative QCD calculations predict the fraction of three jet events i?a(2/cut) as a 
function of y cnt and as- The NNLO QCD prediction for the three-jet rate can be written 
as 

RsiVcut) = a3 - iet( - Vcnt) = a s A 3 (y cut ) + fig B 3 (y cut ) + fi| C 3 (y cnt ), (2) 

Ctat 
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with as = «s(A i )/(2vr) the only free parameter. For this analysis the coefficients A%, B3 
and C3 are taken from [15]. Equation |2] is shown for renormalisation scale \x = Q, where 
Q is the physical scale usually identified with the centre-of-mass energy yfs for hadron 
production in e + e~-annihilation. The terms generated by variation of the renormalisation 
scale parameter = fi/Q are implemented according to [TB] . 

It is well known that for small values of y cnt the fixed order perturbative prediction 
is not reliable, because the expansion parameter (as)L 2 , where L = — \ny cnt , logarithmi- 
cally enhances the higher-order corrections. For instance, (as) ln 2 (0.01) ~ 0(1). Thus, 
one has to perform the all-order resummation of the leading and NLLA contributions. 
This resummation is possible for the Durham algorithm using the coherent branching 
formalism [12l[T6]. Matched NLO+NLLA predictions for jet rates have been compared 
to ALEPH data and good agreement was found, see Fig. 2 of [TTJ . There an improved 
resummation formula was used, where part of the subleading logarithms, that can be 
controlled systematically by including the 'i^-term' [T81420] . are also taken into account 
(NNLO+NLLA+K). In this paper we use the same resummation formula and extend the 
matching to the NNLO accuracy, which requires the expansion of the improved resumma- 
tion formula up to 0(dg). In the resummed prediction we use the one-loop formula for the 
running coupling. One could also use higher-loop running, but the difference would be in 
the coefficients of the subleading logarithms (NNLL and higher), which are not controlled 
systematically in the resummed prediction and therefore are neglected. Expanding the 
resummation formula in [12], we find 

3 2i 

R 3 = Y, ^(CjraaYB^lJ + 0(aj) , (3) 

i=l j=i 

with the i?3 coefficients are given in Table [TJ The dependence on the renormalisation 
scale n can be obtained by making the substitution 

«s -»■ «s ^2 ( ^°" S log n ) ' ( 4 ) 

n=0 ^ 

and keeping the terms up to 0(d|) in the expansion. In Eq. (jlj), 

A) = y c A - 2 -m (5) 

with C*a and C-p are the quadratic Casimir operators of the gauge group in the adjoint 
and fundamental representation, while rif is the number of light flavours (we use ri{—5). 
The matched NNLO+NLLA+K predictions compared to NNLO, NLO and matched 
NLO+NLLA+K, NNLO+NLLA predictions are shown in Fig. Q] 

As opposed to event-shapes, such as the 2/23-distribution [2T| [22]. jet rates do not obey 
simple exponentiation (except for the two-jet rate). For an observable that does not 
exponentiate, the viable matching scheme is the so-called i?-matching [23] . To obtain the 
.R-matched predictions, we subtract the expansion of -R3 from the resummation formula 
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and add the corresponding NNLO prediction given by Eq.2, 

i^ matCh (ycut) = ^(Z/cut) + « S (MV«*) ~ ^ LL (2/cut)) (6) 

+&l{B,{y cut ) - 5 3 NLL (2/cu t )) + &l(c 3 (y cut ) - C 3 NLL (|/ cut )) , 

where A^ LL = C F £j =1 R$' j) V, E 3 NLL = C| £* =2 Rp j ^L j , C 3 NLL = Cf £? =3 i^L*' and 
v4 3 , 5 3 , C 3 as in Eq. |2j Also in the case of jet rates the resummed logarithm is fixed to 
L = — In y cnt unambiguously, and the theoretical uncertainty of the prediction is estimated 
by varying the renormalisation scale (see Fig. [T]) as in [21]. The theoretical uncertainty 
estimated by varying the renormalisation scale is found to be larger by using matched 
NNLO+NLLA+K predictions instead of NNLO+NLLA predictions. 
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Table 1: Coefficients R%'^ of the expansion in Eq. (|2D with .K" = — x — ^?/f, 



3 Analysis Procedure 
3.1 The JADE Detector 

A brief description of the JADE detector focusing on the major detector parts used in 
this analysis is given here. Primarily the momenta of charged and neutral particles are 
used in this analysis. The trajectories of charged particles are mainly reconstructed with 
the central tracking detector, consisting mainly of a large scale jet chamber. The tracking 
detector is located within a 0.48 T solenoidal magnetic field. Charged and neutral particles 
except muons and neutrinos are reconstructed with an electromagnetic calorimeter which 
surrounds the magnetic coil. The calorimeter consists of lead glass Cherenkov counters 
and is separated in a barrel and two endcap regions. A more complete summary of the 
JADE detector can be found in |25j . 
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Fig. 1: The QCD predictions are shown for the Durham three-jet rate calculated with 
a value of as(Mzo)=0.1180 and at a centre-of-mass energy of 35 GeV. In a) the NNLO 
prediction is compared to the matched NNLO+NLLA and NNLO+NLLA+K predictions. 
The insert shows the ratio between the different predictions. In b) the matched and 
unmatched NLO predictions are compared to the matched NNLO prediction with the 
insert presenting the ratio to the fixed order and the matched NLO predictions. For all 
calculations the uncertainty band reflects the uncertainty originating from setting the 
renormalisation scale factor = /x/Q to 0.5 and 2. 



3.2 Data and Monte Carlo Samples 

For this analysis we use the identical data sample as used for previous JADE analyses [HI 
|2"4"1T26] . The data were taken between 1979 and 1986, adding up to an integrated luminosity 
of about 195 pb" 1 . It is subdivided into six data taking periods with different average 
centre-of-mass energies of 14 GeV, 22 GeV, 34.6 GeV, 35 GeV, 38.3 GeV and 43.8 GeV. 
The number of selected events ranges between 1403 events at 22 GeV and 20876 events 
at 35 GeV. 

To correct for acceptance and resolution effects as well as for hadronisation effects a 
large sample of Monte Carlo events is generated. The Monte Carlo generators are tuned to 
match events taken with the OPAL experiment at a centre-of-mass energy of 91 GeV [271 
128] . Using these same parameters, except of setting the appropriate centre-of-mass energy, 
a good description of the JADE data, down to smallest energies, is achieved [29J. For the 
correction of acceptance and resolution effects the events are passed through a simulation 
of the JADE experiment. Events generated with PYTHIA 5.7 [30J are used as default 
for the correction of acceptance and resolution effects. As a cross check events simulated 
with the HERWIG 5.9 [31] event generator are utilised. To assess the changes in the 
three-jet rate distribution originating from the transition from partons to hadrons three 
different Monte Carlo generators are applied. As default PYTHIA 6.158 is used and 
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events produced with the HERWIG 6.2 and ARIADNE 4.11 [32] generators are used as 
a consistency check. 

Due to the good agreement between data and the predictions from the PYTHIA, 
HERWIG and ARIADNE Monte Carlo event generators (see section loTTj) we do not utilise 
event generators which incorporate high-multiplicity matrix elements which were not yet 
tuned to our purpose. This choice is justified by the fact that generators used in our study 
use leading-order matrix elements combined with a leading-logarithm parton shower which 
do provide a satisfactory description of the three parton final states studied in this analysis 
as shown in [2"HI3"3"]. 

3.3 Event Selection 

The measurement of the strong coupling as is based on the analysis of well measured 
hadronic events. A detailed description of the hadronic event selection can be found 
in [H]. Cuts applied to the events are based on the number of charged tracks, total 
visible energy and momentum imbalance. The cut on the momentum imbalance reduces 
events emitting a high energetic photon in the initial state, leading to a reduced hadronic 
centre-of-mass energy. The cuts on the number of charged tracks and total visible energy 
minimise to an insignificant level the number of events from hadronic tau decays and from 
hadronic final states originating from two-photon scattering. 

3.4 Corrections to the Data 

The corrections applied to the data follow exactly the same procedure as summarised 
in [14] . For the calculation of the three-jet rate all charged tracks and electromagnetic 
clusters are considered. The estimated minimum ionizing energy from tracks associated 
with electromagnetic calorimeter clusters is subtracted from the cluster energies. For the 
correction procedure two different categories of jet rate distributions are defined for sim- 
ulated events, the so-called detector level and the hadron level distribution. The detector 
level distribution of simulated events is obtained by using all selected tracks of charged 
particles and electromagnetic clusters. The hadron level distribution is obtained by using 
the true four-momenta of stable particles, where particles with a lifetime of r > 300 ps are 
declared as stable. Simulated events with photon initial state radiation (ISR) leading to 
the centre-of-mass energy reduced by more than 0.15 GeV are rejected from the hadron 
leveVQ. Thus the correction for experimental effects explained below also takes care of 
residual effects due to initial state radiation. Only hadronic events originating from the 
primary production of u,d,s or c quarks are considered. 

Before correcting the three-jet rate data distribution the expected contribution from 
e + e~ — > bb-events as expected from simulated events at the detector level is subtracted 
from the three jet-rate. About 1/11 of all qq-events are bb-events and the expected num- 
ber of bb-events is subtracted from the observed number of data events at each y cut -bin. 
The hadronic events used in this analysis correspond to events with e + e~ annihilating 
to a pair of u,d,s or c-quarks. The distribution corrected for e + e~ — > bb-events is then 

2 The cut value of 0.15 GeV is purely technical and corresponds to a clean separation of the soft and 
hard ISR events in the simulation. 
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multiplied bin-by-bin with the ratio of the hadron level distribution divided by detector 
level distribution. A correction method based on a matrix unfolding method returns com- 
patible results within statistical uncertainties [29J. The impact on the measurement due 
to changes of the 6-quark fragmentation in the simulation are covered by the systematic 
uncertainty (see section H]) assigned to the correction for bbevents • The numerical 
results of the corrected distributions are summarised in Tables [6] and [71 

4 Systematic Uncertainties 

In order to assess the systematic uncertainty of the measurement of the strong coupling 
as we evaluate several possible sources. For each variation the difference to the result 
with respect to the default analysis procedure is taken as a contribution to the systematic 
uncertainty. The systematic uncertainty is assumed to be symmetric around the default 
value. No systematic uncertainty is evaluated related to the fact that massless theoretical 
predictions are used since the contribution from e + e~ — > bb-events is subtracted from the 
data distribution (see section [33]). The uncertainty originating from correcting for ISR 
effects is small and no systematic uncertainty is assigned. 

4.1 Experimental Uncertainties 

The assessment of the experimental uncertainty follows exactly the procedure described in 
detail in [14J . The analysis is repeated with a slightly modified event and track selection, 
using a different reconstruction software, using different MC models for the correction of 
detector effects and modified fit ranges. In addition the amount of subtracted e + e~ — > bb- 
events is modified by ±5%. The differences between the results obtained from the modified 
fits and the default fits are added in quadrature and taken as the combined experimental 
systematic uncertainty. The main contribution originates from using a different version 
of the reconstruction software and using HERWIG instead of PYTHIA for the correction 
of detector effects. 

4.2 Hadronisation 

The standard analysis uses PYTHIA to evaluate the change of the three-jet rate distribu- 
tion originating from the transition from partons to stable hadrons (see section I5T2|) . Only 
the process e + e~ to a pair of u,d,s or c-quarks is simulated. To assess the uncertainty 
the fit is repeated with alternative Monte Carlo models. For this we use HERWIG and 
ARIADNE instead of PYTHIA and the difference to the standard PYTHIA correction 
is taken as systematic uncertainty. In all cases the larger difference is seen by using the 
HERWIG Monte Carlo simulation instead of PYTHIA. 

A variation of the parameters describing the hadronization model leads to significant 
smaller changes on the measurement of the strong coupling as than applying a different 
hadronization scheme, like HERWIG [29J. For this reason we quote only the largest 
uncertainty coming from using a different hadronization model. As described in 13.21 
we use different versions of the Monte Carlo generator for detector and hadronisation 
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corrections. However, the hadronisation corrections for each generator are consistent 
within their statistical uncertainty. 

4.3 Theoretical Uncertainties 

The theoretical prediction of the three-jet rate distribution is a truncated asymptotic 
power series. The uncertainties originating from missing higher order terms are assessed 
by varying the renormalisation scale parameter x ii = \ij \fs. For this x M is set to two and to 
0.5. The larger deviation from the fit using the default setting x M = 1 is taken as system- 
atic uncertainty. We assess the theoretical uncertainties by applying NNLO+NLLA+K 
QCD predictions in the fit, since the uncertainties using NNLO+NLLA predictions with- 
out fT-term are found to be smaller. Effects from electroweak corrections are neglected. 



5 Results 

5.1 Three- Jet Rate Distributions 

The three-jet rates as a function of y cut at y/s= 14, 22, 34.6, 35, 38.3 and 43.8 GeV 
are shown in Fig. [2j The rates are corrected for resolution and acceptance effects. The 
estimated contributions from e + e~ — > bb-events are subtracted. The distributions are 
compared with the predictions obtained from Monte Carlo models used in this analysis. 
All models reproduce the data distribution well. Most of the data points are within the 
one sigma uncertainty band and a couple of points show a deviation of up to two sigma 
uncertainty. There is an apparent deviation of the simulation from the data in Fig. [2] 
(bottom left) for i/s=38.3 GeV. This is due to the positive correlations between the data 
points which implies that a fluctuation in one data point is also visible in the neighboring 
data points. For additional information inserts show the deviation of the three-jet rate 
obtained from Monte Carlo with respect to the data, normalised to the statistical and 
experimental uncertainty added in quadrature!!. 



3 Please note that correlations between the points are present and not taken into account for the insert 
plot. 
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Fig. 2: The three-jet rate as a function of the resolution parameter y cut is shown for 
centre-of-mass energies y^s of 14.0 to 43.8 GeV. The data distributions are corrected for 
resolution and acceptance effects from the detector and the contributions from e + e _ — > bb- 
events are subtracted. Data points are shown with statistical uncertainty (inner part) and 
combined statistical and experimental uncertainty. The histograms show the comparison 
to the predictions obtained with Monte Carlos simulation using PYTHIA, HERWIG and 
ARIADNE. The inserts show the deviation from the simulated distribution normalised to 
the combined statistical and experimental uncertainty. 
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5.2 Determination of cxs 

The value of the strong coupling is determined by a minimum-^ 2 fit of the matched 
NNLO+NLLA+K predictions to the corrected data distributions, separately for each 
centre-of-mass energy. The reduced renormalisation scale parameter x M = \xj\fs in the 
theoretical predictions is set to the natural choice The QCD predictions describe 

the three-jet rate at parton level only. To correct for hadronisation effects the matched 
QCD predictions are multiplied at each y cut point by the ratio of the hadron level dis- 
tribution divided by the parton level distribution obtained from simulated events. The 
parton level distribution in simulation is obtained from the final state partons after the 
parton shower has terminated, i.e. just before the hadronisation step. The hadron level 
is defined in an identical way as described in section I3.4[ containing only hadronic events 
with e + e~ annihilating to a pair of u,d,s or c-quarks. The ratio between the hadron level 
divided by the parton level estimated with simulated events is shown in the appendix in 
Fig. [HI Using the QCD-prediction corrected for hadronisation effects and the corrected 
data distribution a x 2 -value is calculated 

n 

X 2 = X> 3li - Riasf^WiRsr'URs, - R(a s )tr) ( 7 ) 

i,3 

with i and j being the y cut points in the chosen fit range and the i?(as)3 h f° vames are the 
predicted values of the three-jet rate. Each event can contribute to several points in the 
three-jet rate distribution leading to correlation between different y cut points. For this 
reason the covariance matrix is not diagonal and the off-diagonal elements have to be 
computed. We follow the approach described in [14J using 1000 subsamples with 1000 
simulated events each. 

When choosing the fit range, we took the following considerations into account. The 
corrections applied to the data (see section I3.4|) reverting the imperfectness of detector 
and the correction of the Q CD-predictions due to hadronisation effects are required to 
be small. In addition we require the leading log contribution in the low y cut region of 
the three-jet rate distribution to be well below unity to ensure that the NLLA is valid. 
The leading log term is proportional to as • ln 2 (l/y cu t) and requiring this term to be 
well below unity leads us to a lower limit of y cu t= 0.01 assuming a value of the strong 
coupling ag(M z o)=0.118. The upper limit is determined by requiring the leading order 
contribution A(y cut ) to be larger than zero. The corrections are small in this range with 
the detector corrections being less than 30% and the hadronisation corrections being less 
than 30%, apart from the corrections at the centre-of-mass energy of 14 GeV, which are 
up to 70%. The considerations described above lead to a fit range from 0.01 to 0.2, which 
is identical to the fit range used for the determination of as using the differential 2/23 event 
shape distribution [2~4] . 
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Fig. 3: The fit results at centre-of-mass energies from 14.0 to 43.8 GeV are shown. The 
inserts show the deviation of the data points from the QCD-prediction with the a s -value 
obtained from the fit, normalised to the combined statistical and experimental error. 
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The results of the fits to the three-jet rate distribution at the various centre-of-mass 
energies are shown in Fig. |3j The numerical results of the fits are summarised in Table EJ 
The statistical uncertainty is obtained from the fit and the systematic uncertainty is 
evaluated as described in section HI 

The fitted theory generally describes the data at hadron level well within the fit ranges, 
as seen in Fig. [3] and confirmed by the x 2 /d.o.f. values in Table EJ The x 2 /d-o.f. values 
are based only on the statistical uncertainties of the data and thus it is reasonable that 
they are larger than unity. The extrapolations outside of the fit regions also provide a 
reasonable description of the data. 

For comparison the numerical results for fits using NLO+NLLA+K predictions are 
compiled in Table [3J The values of as are consistently smaller and the theoretical un- 
certainty is about a factor of three larger compared to the fit using NNLO+NLLA+K 
calculations. The fit quality reflected by the x 2 /d.o.f .-value is similar for both QCD pre- 
dictions, but slightly better for the fit based on NNLO+NLLA+K predictions. The differ- 
ence between the as value returned with matched NLO and matched NNLO predictions 
is of similar size as the theoretical uncertainty evaluated for the fit using NLO+NLLA+K 
calculations. The reduction of the theoretical uncertainty using NNLO+NLLA+K QCD 
predictions is associated to the higher-order terms available in the new calculations. 

The result obtained here cannot be compared directly to the measurement of as us- 
ing the identical data set and the y 2 3 event shape observable [24]. While for the as 
measurement in [23] matched NNLO+NLLA QCD predictions are used, in this analysis 
NNLO+NLLA+K predictions are applied, which take subleading logarithms into account. 



[GeV] 




stat. 


exp. 


hadr. 


scale 


X 2 /d.o.f. 


14.00 


0.1704 


0.0029 


0.0019 


0.0079 


0.0028 


62.95/10 


22.00 


0.1562 


0.0044 


0.0030 


0.0128 


0.0013 


13.63/10 


34.60 


0.1399 


0.0015 


0.0023 


0.0086 


0.0009 


13.17/10 


35.00 


0.1469 


0.0013 


0.0043 


0.0086 


0.0011 


11.69/10 


38.30 


0.1375 


0.0052 


0.0084 


0.0096 


0.0013 


37.58/10 


43.80 


0.1329 


0.0029 


0.0032 


0.0045 


0.0008 


20.29/10 



Table 2: The value of as using matched NNLO+NLLA+K predictions together with 
the statistical, experimental, hadronisation and theoretical uncertainties as described in 
section HI The last column shows the xVd-o.f. value of the fit obtained with statistical 
uncertainties only at the respective energy points. 



5.3 Combination of as Measurements 

The results of the measurements of the strong coupling as at the various centre-of-mass 
energies are combined to a single value of as(M z o). For this the values of cts(\/s) obtained 
are evolved to a common energy scale Mgo and combined using a weighted mean. The 
theoretical uncertainty as well as the uncertainty originating from modelling the hadroni- 
sation process are likewise determined by calculating the weighted mean of the uncertainty 
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[GeV] 




stat. 


exp. 


hadr. 


scale 


X 2 /d.o.f. 


14.00 


0.1648 


0.0026 


0.0020 


0.0068 


0.0039 


64.44/10 


22.00 


0.1518 


0.0040 


0.0026 


0.0118 


0.0043 


13.12/10 


34.60 


0.1367 


0.0014 


0.0022 


0.0082 


0.0025 


12.39/10 


35.00 


0.1432 


0.0012 


0.0039 


0.0083 


0.0028 


14.54/10 


38.30 


0.1338 


0.0048 


0.0075 


0.0089 


0.0042 


38.47/10 


43.80 


0.1300 


0.0027 


0.0029 


0.0043 


0.0021 


21.65/10 



Table 3: The value of as using NLO+NLLA+K predictions together with the statistical, 
experimental, hadronisation and theoretical uncertainty as described in section [4] . The 
last column shows the x 2 /d.o.f. value of the fit obtained with statistical uncertainties only 
at the respective energy points. 

evaluated at each single energy point. The difficulties arise in estimating the correlations 
between the systematic uncertainties obtained at the different energy points. The identi- 
cal problem is present for the combination of as results from the LEP-collaborations and 
for this reason we use the same method as outlined in [T3JI2UI31]- For this analysis we 
combine the values obtained at centre-of-mass energies between 22 and 43.8 GeV, exclud- 
ing the result obtained at 14 GeV because of the large hadronisation corrections and the 
large xVd.o.f. returned by the fit. The combination of the results obtained between 22 
and 43.8 GeV results in 

a s (M z o) = 0.1199 ± 0.0010(stat.) ± 0.0021(exp.) ± 0.0054(had.) ± 0.0007(theo.) 

consistent with the world average of as = 0.1184 ± 0.0007 [35]. To the combined 
result the value of as measured at 22 GeV contributes with a weight of 0.13, at 34.6 GeV 
with 0.23, at 35.0 GeV with 0.16, at 38 GeV with 0.07 and at 43.8 GeV with a weight 
of 0.41. The results at the various centre-of-mass energies are visualised in Fig. HJ Here, 
the values obtained at 34.6 and 35.0 GeV are combined to a single value using the same 
method as described above 

a s (34.8 GeV) = 0.1431 ± 0.0010(stat.) ± 0.0028(exp.) ± 0.0086(had.) ± 0.0010(theo.). 
The combination of the as measurements at centre-of-mass energies between 22 and 
43.8 GeV determined with NLO+NLLA+K predictions leads to 

a s (M z o) = 0.1175 ± 0.0010(stat.) ± 0.0020(exp.) ± 0.0052(had.) ± 0.0020(theo.). 

5.4 Simultaneous Variation of a§ and the Renormalisation Scale 

Besides fixing the renormalisation scale parameter to the natural choice x M =l we repeat 
the fit as a cross-check with both, as and x M , being varied within the minimisation pro- 
cedure, thus choosing the so called optimal renormalisation scheme |36j. The results of 
the fits are summarised in Table HI The x 2_vames obtained are almost identical to that 
from the default fit. The results for the scale x M =a;° pt at the minimal xVd.o.f. -value vary 
between 0.31 and 2.40, all being within the errors consistent with the x M -range used in 
the systematic variation of the default fit. To estimate the theoretical uncertainty the fit 
is repeated with being set to twice and half of the value obtained for the optimal scale 
x° pt . The as-results at the various energy points between 22 and 43.8 GeV are combined 
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Fig. 4: The open and solid points show the measurements of a$ from the three-jet rate 
at the various centre-of-mass energies. The error bars indicate the statistical (inner part) 
and the total error. The lines show the world average of as(M z o) [35J. The results at 34.6 
and 35.0 GeV are combined to a single value. The result obtained at 14 GeV (open point) 
is not used as input for the combined value of as- The result shown at a centre-of-mass 
energy of 91.2 GeV (triangle) is obtained from a fit to the three-jet rate using NNLO 
predictions only with data taken by the ALEPH detector |13j . 



to a single value using the method described in 15.31 The combined value is: 

a s (M z o) = 0.1204 ± 0.0009(stat.) ± 0.0021(exp.) ± 0.0059(had.) ± 0.0008(theo.), 
leading to almost the identical result as obtained with the default fit. The variation 
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of as an d x 2 /d-o.f. with respect to the renormalisation scale is summarised in Fig. |5j 
The changes for both, the strong coupling as as well as the x 2 /d.o.f. -value, are small 
within the a^-range considered, indicating that our results are only moderately sensitive 
to missing higher order terms. 



[GeV] 


(*s(y/s) 


stat. 


exp. 


hadr. 


scale 


™opt 
X M 


Corr. 


X 2 /d.o.f. 


14.00 


0.1739 


0.0032 


0.0020 


0.0088 


0.0149 


0.32 ±0.06 


0.26 


57.04/9 


22.00 


0.1559 


0.0044 


0.0029 


0.0127 


0.0014 


1.11 ± 1.07 


-0.44 


13.62/9 


34.60 


0.1390 


0.0015 


0.0023 


0.0086 


0.0007 


2.40 ± 2.08 


-0.39 


12.89/9 


35.00 


0.1478 


0.0014 


0.0044 


0.0088 


0.0012 


0.58 ±0.29 


-0.49 


11.03/9 


38.30 


0.1387 


0.0054 


0.0087 


0.0100 


0.0016 


0.55 ±0.36 


0.15 


37.10/9 


43.80 


0.1337 


0.0029 


0.0034 


0.0048 


0.0042 


0.31 ±0.16 


0.27 


18.82/9 



Table 4: The value of as using matched NNLO±NLLA±K predictions applying the opti- 
mised renormalisation scheme. The result is shown together with the statistical, experi- 
mental, hadronisation, theoretical uncertainties, the renormalisation scale parameter x° pt , 
the correlation between as and x° pt as well as the x 2 /d.o.f. value of the fit. 
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Fig. 5: The variation of the strong coupling as as a function of the renormalisation scale 
parameter and the corresponding x 2 /d.o.f .-value for the data taken at a centre-of- 
mass energies from 14 to 43.8 GeV using the matched NNLO+NLLA+K predictions. 
The arrows indicate the variation of the renormalisation scale parameter used to evaluate 
the theoretical systematic uncertainty. 
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5.5 Measurements of a§ using NNLO Predictions only 

To compare our result with results obtained with NNLO predictions only, we repeat the 
fit without NLLA matching. The fit range is chosen to be identical to the fit range used in 
the default fit. We perform two different fits using different choices for the renormalisation 
scale parameter: a fixed value of x M =l and the optimised renormalisation scheme, where 
x M is an additional free parameter in the fit. The numbers in Table [5] show the result 
obtained at a centre-of-mass energy of 35.0 GeV, the energy point with the largest number 
of selected hadronic events. The results obtained at all energy points are summarised in 
the appendix in Table M and Table M Fig- M shows the fit result using NNLO predictions 
with and x A1 =x° pt at a centre-of-mass energy of 35.0 GeV. 

The shape of the three-jet distribution is not well matched by using NNLO predictions 
only, even at large y cut values. The x 2 /d.o.f. -value obtained with x M =l increases signifi- 
cantly compared to the corresponding NNLO+NLLA+K fit. The fit using the optimised 
renormalisation returns reasonable x 2 /d.o.f. -values and smaller x M -values. A similar be- 
havior was already observed using NLO predictions only with the renormalisation factor 
being set to x M =l and £ M =a;° pt [SUEZ] wnere variation of the renormalisation scale fac- 
tor led to an improved description of the data as well. For both renormalisation scale 
schemes the scale uncertainty is considerably increased compared to the measurement 
using matched NNLO+NLLA+K predictions. In addition the difference between the fit 
using the natural and the optimised renormalisation scale is increased compared to a com- 
parison using matched NNLO+NLLA+K predictions. A large sensitivity to the choice of 
the fit range is observed for a fit using NNLO predictions only. Again, the measurement 
of as with data taken at a centre-of-mass energy of 14 GeV returns by far the largest 

rVd.cf.. 

We conclude, contrary to [13] which performs a fit to a single y cnt -bm only and there- 
fore being insensitive to the shape of the distribution, that resummation affects the fit 
significantly by decreasing the scale dependence and making the result of the fit more 
reliable. For this reason we consider the result based only on NNLO predictions as a 
cross-check only. 





a s (35 GeV) 


stat. 


exp. 


hadr. 


scale 




X 2 /d.o.f. 


rp ™Opt 


0.1426 
0.1486 


0.0012 
0.0014 


0.0040 
0.0049 


0.0089 
0.0089 


0.0015 
0.0050 


1.0 

0.23 + 0.03 


45.30/10 
11.42/9 



Table 5: The value of as determined with NNLO predictions only together with the 
statistical, experimental, hadronisation and theoretical uncertainty as described in section 
HJ The last columns show the value for x M and the xVd.o.f. The first row shows the fit 
result for a fixed renormalisation scale parameter x M =l, the second row for x M being varied 
within the fit. The x 2 - value is obtained with the statistical uncertainty only, taking bin- 
to-bin correlations into account. 
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Fig. 6: The three-jet rate distribution at a centre-of-mass energy of 35 GeV is shown 
together with the NNLO QCD prediction evaluated at the as-value obtained from the 
least-x 2 fit (solid line). In addition the QCD-predictions obtained with the optimised 
renormalisation scale parameter (dash-dotted line) is shown. The QCD predictions are 
corrected for hadronisation effects. The insert shows the deviation of the data points from 
the QCD-prediction with the a§- value obtained from the fit, normalised to the combined 
statistical and experimental error. 

6 Summary 

In this paper we present a measurement of the strong coupling a s using the three- 
jet rate taken with the JADE experiment at centre-of-mass energies between 14 and 
43.8 GeV. The three-jet rate is compared to simulated events obtained with PYTHIA, 
HERWIG and ARIADNE. All Monte Carlo models reproduce the measured distribu- 
tions well. The strong coupling as is measured by a fit to matched NNLO+NLLA+K 
predictions. The combined value using the measurements between 22 and 43.8 GeV 
results in as(M z o) = 0.1199 ± 0.0060 (total error). The value is consistent with the 
world average [35]. The theory uncertainty has shrunk considerably by using matched 
NNLO+NLLA+K calculations and among the uncertainties it is the smallest. The dom- 
inant uncertainty originates from applying different Monte Carlo models to estimate the 
transition from partons to hadrons. A fit using NNLO predictions only with the renor- 
malisation scale parameter x M set to one cannot describe the shape of the three-jet rate 
distribution well. The shape can only be described well if as and are fitted simultane- 
ously. 

Fig. [7] shows a comparison of the as measurement presented in this paper with previous 
as measurements using higher QCD predictions. The results obtained with a fit to the 
l/23-distribution using data taken with the JADE [24] and the OPAL experiments [33] are 
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Fig. 7: Comparison of the «s measurement obtained in this analysis with previous as 
measurements using the three-jet rate, the four-jet rate and the 2/23-differential distribu- 
tion in e + e~-annihilation between 14 and 91 GeV using matched and unmatched QCD 
predictions. The 1 st measurement (top down) is the result of the present analysis, the 2 nd 
is the measurement using the differential two-jet rate distribution 2/23 with NNLO+NLLA 
calculations and data taken between 14 and 44 GeV [23], the 3 rd is the measurement 
using the four-jet rate using matched NLO+NLLA+K predictions and data between 14 
and 44 GeV [Hj, the 4 th is the measurement using the differential two-jet rate distribu- 
tion with NNLO+NLLA calculations and data taken at 91 GeV [33j and the 5 th is a 
measurement using the three-jet rate at a centre-of-mass energy of 91 GeV using NNLO 
predictions only [T5] . 

expected to be strongly correlated with the three-jet rate and are considered as a good 
cross check. The measured mean values, similar to the result obtained with this analysis, 
are slightly above the world average value. The measurements using the four-jet rate and 
data taken with the JADE experiment [H] and the three-jet rate using data taken with 
the ALEPH experiment [13] are consistent within the uncertainties. All measurements 
are in agreement with the world average value of «s [35] . 
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logio(ycut) 


i? 3 (14 GeV) 


#3(22 GeV) 


#3(34.6 GeV) 


-4.93 








-4.80 








-4.68 








-4.55 








-4.43 






0.002 ±0.000 ±0.003 


-4.30 






0.001 ±0.000 ±0.001 


-4.18 






0.001 ±0.000 ±0.001 


-4.05 




0.002 ±0.001 ±0.002 


0.002 ±0.000 ±0.003 


-3.93 




0.001 ±0.001 ±0.001 


0.002 ±0.000 ±0.001 


-3.81 




0.001 ±0.001 ±0.001 


0.002 ±0.000 ±0.001 


-3.68 




0.000 ± 0.000 ± 0.004 


0.002 ±0.000 ±0.001 


-3.56 


0.001 ±0.001 ±0.001 


0.000 ±0.001 ±0.000 


0.003 ±0.001 ±0.001 


-3.43 


0.002 ±0.001 ±0.002 


0.001 ±0.001 ±0.001 


0.006 ±0.001 ±0.002 


-3.31 


0.002 ±0.001 ±0.001 


0.003 ±0.002 ±0.002 


0.014 ±0.001 ±0.002 


-3.18 


0.003 ±0.001 ±0.001 


0.006 ±0.002 ±0.002 


0.026 ±0.001 ±0.002 


-3.06 


0.005 ±0.002 ±0.001 


0.013 ±0.003 ±0.002 


0.053 ±0.002 ±0.005 


-2.93 


0.008 ±0.002 ±0.001 


0.025 ±0.004 ±0.003 


0.098 ±0.003 ±0.006 


-2.81 


0.015 ±0.003 ±0.003 


0.047 ±0.006 ±0.003 


0.167 ±0.003 ±0.007 


-2.68 


0.027 ±0.004 ±0.007 


0.103 ±0.008 ±0.009 


0.251 ±0.004 ±0.011 


-2.56 


0.061 ±0.005 ±0.009 


0.168 ±0.010 ±0.013 


0.331 ±0.004 ±0.012 


-2.43 


0.109 ±0.007 ±0.014 


0.252 ±0.012 ±0.018 


0.393 ±0.004 ±0.013 


-2.31 


0.188 ±0.009 ±0.014 


0.343 ±0.013 ±0.026 


0.414 ±0.004 ±0.009 


-2.18 


0.281 ±0.010 ±0.009 


0.424 ±0.014 ±0.015 


0.401 ±0.004 ±0.005 


-2.06 


0.374 ±0.011 ±0.019 


0.436 ±0.014 ±0.022 


0.369 ±0.004 ±0.005 


-1.93 


0.461 ±0.011 ±0.022 


0.413 ±0.014 ±0.010 


0.329 ±0.004 ±0.009 


-1.80 


0.485 ±0.011 ±0.015 


0.343 ±0.013 ±0.017 


0.286 ±0.004 ±0.005 


-1.68 


0.448 ±0.011 ±0.007 


0.303 ±0.013 ±0.015 


0.248 ±0.004 ±0.007 


-1.55 


0.360 ±0.011 ±0.010 


0.240 ±0.012 ±0.010 


0.207 ±0.004 ±0.008 


-1.43 


0.281 ±0.010 ±0.011 


0.190 ±0.011 ±0.010 


0.171 ±0.003 ±0.006 


-1.30 


0.205 ±0.009 ±0.013 


0.154 ±0.010 ±0.012 


0.132 ±0.003 ±0.006 


-1.18 


0.132 ± 0.008 ± 0.016 


0.115 ± 0.009 ± 0.008 


0.100 ± 0.003 ± 0.005 


-1.05 


0.084 ±0.006 ±0.012 


0.078 ±0.007 ±0.006 


0.068 ±0.002 ±0.003 


-0.93 


0.045 ±0.005 ±0.011 


0.052 ±0.006 ±0.008 


0.044 ±0.002 ±0.001 


-0.81 


0.012 ±0.002 ±0.007 


0.030 ±0.005 ±0.006 


0.024 ±0.001 ±0.002 


-0.68 


0.005 ±0.002 ±0.002 


0.014 ±0.003 ±0.002 


0.009 ±0.001 ±0.001 


-0.56 


0.003 ±0.001 ±0.002 




0.001 ±0.000 ±0.000 



Table 6: Hadron-level value of the three-jet fraction using the Durham algorithm at 14, 
22 and 34.6 GeV. The value is corrected for contributions from e + e~ — > bb-events. In all 
cases the first quoted error indicates the statistical error while the second quoted error 
corresponds to the total experimental uncertainty. Uncertainties consistent with zero 
indicate that the corresponding value is smaller than the precision shown in the Table. 
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logio(ycut) 


i? 3 (35 GeV) 


#3(38.3 GeV) 


i? 3 (43.8 GeV) 


-4.93 








-4.80 








-4.68 








-4.55 








-4.43 








-4.30 








-4.18 






0.003 ±0.001 ±0.006 


-4.05 






0.005 ±0.001 ±0.006 


-3.93 


0.001 ±0.000 ±0.001 


0.001 ±0.001 ±0.000 


0.003 ±0.001 ±0.003 


-3.81 


0.001 ±0.000 ±0.001 


0.001 ±0.001 ±0.001 


0.003 ±0.001 ±0.001 


-3.68 


0.002 ±0.000 ±0.001 


0.003 ±0.002 ±0.002 


0.005 ±0.001 ±0.001 


-3.56 


0.003 ±0.000 ±0.001 


0.005 ±0.002 ±0.002 


0.007 ±0.001 ±0.001 


-3.43 


0.006 ±0.001 ±0.001 


0.008 ±0.003 ±0.003 


0.012 ±0.002 ±0.002 


-3.31 


0.011 ±0.001 ±0.001 


0.018 ±0.004 ±0.003 


0.026 ±0.003 ±0.003 


-3.18 


0.024 ±0.001 ±0.001 


0.040 ± 0.005 ± 0.007 


0.053 ±0.004 ±0.004 


-3.06 


0.050 ±0.002 ±0.003 


0.064 ±0.007 ±0.005 


0.101 ±0.005 ±0.007 


-2.93 


0.096 ±0.002 ±0.002 


0.123 ±0.009 ±0.009 


0.172 ±0.006 ±0.007 


-2.81 


0.165 ±0.003 ±0.004 


0.187 ±0.011 ±0.010 


0.254 ±0.007 ±0.011 


-2.68 


0.245 ± 0.003 ± 0.008 


0.270 ±0.012 ±0.014 


0.339 ±0.008 ±0.006 


-2.56 


0.326 ±0.004 ±0.009 


0.346 ±0.013 ±0.009 


0.390 ±0.008 ±0.008 


-2.43 


0.381 ±0.004 ±0.003 


0.379 ±0.013 ±0.006 


0.415 ±0.008 ±0.009 


-2.31 


0.401 ±0.004 ±0.004 


0.386 ±0.013 ±0.012 


0.411 ±0.008 ±0.007 


-2.18 


0.397 ±0.004 ±0.004 


0.365 ±0.013 ±0.010 


0.381 ±0.008 ±0.007 


-2.06 


0.377 ±0.004 ±0.004 


0.328 ±0.013 ±0.014 


0.346 ±0.008 ±0.009 


-1.93 


0.343 ±0.004 ±0.008 


0.300 ±0.013 ±0.013 


0.307 ±0.008 ±0.005 


-1.80 


0.304 ±0.004 ±0.007 


0.279 ±0.012 ±0.015 


0.281 ±0.007 ±0.011 


-1.68 


0.261 ±0.003 ±0.011 


0.268 ±0.012 ±0.017 


0.242 ±0.007 ±0.006 


-1.55 


0.220 ±0.003 ±0.011 


0.239 ±0.012 ±0.021 


0.195 ±0.006 ±0.006 


-1.43 


0.176 ±0.003 ±0.009 


0.198 ±0.011 ±0.017 


0.156 ±0.006 ±0.006 


-1.30 


0.140 ±0.003 ±0.008 


0.152 ±0.010 ±0.015 


0.126 ±0.005 ±0.007 


-1.18 


0.106 ± 0.002 ± 0.006 


0.117 ± 0.009 ± 0.012 


0.097 ± 0.005 ± 0.004 


-1.05 


0.077 ±0.002 ±0.006 


0.088 ±0.008 ±0.009 


0.072 ±0.004 ±0.005 


-0.93 


0.049 ±0.002 ±0.001 


0.060 ±0.006 ±0.010 


0.043 ±0.003 ±0.004 


-0.81 


0.028 ±0.001 ±0.003 


0.034 ±0.005 ±0.007 


0.023 ±0.003 ±0.001 


-0.68 


0.011 ±0.001 ±0.002 


0.007 ±0.002 ±0.002 


0.008 ±0.002 ±0.002 


-0.56 


0.002 ±0.000 ±0.000 


0.006 ±0.002 ±0.009 


0.003 ±0.001 ±0.002 



Table 7: Hadron-level value of the three-jet fraction using the Durham algorithm at 35, 
38.3 and 43.8 GeV. The value is corrected for contributions from e + e~ — > bb-events. In 
all cases the first quoted error indicates the statistical error while the second quoted error 
corresponds to the total experimental uncertainty. Uncertainties consistent with zero 
indicate that the corresponding value is smaller than the precision shown in the Table. 
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B Hadronisation correction estimated with simulated 
events 
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Fig. 8: The ratio of the hadron level prediction divided by the partem level prediction as 
a function of the resolution parameter y cut shown for centre-of-mass energies y/s between 
14 and 43.8 GeV. 
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C Fit results with NNLO predictions 



[GeV] 


as(v^) 


stat. 


exp. 


hadr. 


scale 


X 2 /d.o.f. 


14.00 


0.1629 


0.0023 


0.0036 


0.0047 


0.0028 


80.24/10 


22.00 


0.1499 


0.0037 


0.0039 


0.0105 


0.0018 


18.40/10 


34.60 


0.1364 


0.0014 


0.0029 


0.0087 


0.0013 


32.59/10 


35.00 


0.1426 


0.0012 


0.0040 


0.0089 


0.0015 


45.30/10 


38.30 


0.1305 


0.0044 


0.0064 


0.0085 


0.0015 


47.55/10 


43.80 


0.1301 


0.0026 


0.0031 


0.0050 


0.0010 


25.48/10 



Table 8: The value of «s and the statistical, experimental, hadronisation and theoretical 
uncertainty as described in section H] using NNLO predictions only. The last column shows 
the x 2 /d.o.f. value of the fit for all energy points. 



[GeV] 


as(v^) 


stat. 


exp. 


hadr. 


scale 


™opt 
X W 


Corr. 


X 2 /d.o.f. 


14.00 


0.1748 


0.0033 


0.0031 


0.0094 


0.0174 


0.22 ±0.02 


-0.68 


61.29/9 


22.00 


0.1547 


0.0042 


0.0042 


0.0123 


0.0045 


0.28 ±0.12 


-0.77 


16.01/9 


34.60 


0.1407 


0.0016 


0.0031 


0.0087 


0.0053 


0.23 ±0.05 


-0.75 


16.83/9 


35.00 


0.1486 


0.0014 


0.0049 


0.0088 


0.0050 


0.23 ±0.03 


-0.72 


11.42/9 


38.30 


0.1651 


0.0053 


0.0090 


0.0068 


0.0274 


0.12 ±0.01 


0.58 


23.98/9 


43.80 


0.1396 


0.0037 


0.0050 


0.0039 


0.0068 


0.14 ±0.04 


-0.71 


16.66/9 



Table 9: The value of as and the statistical, experimental, hadronisation, renormalisation 
scale parameter x° pt , the correlation between as and the renormalisation scale parameter 
and the x 2 /d.o.f. value of the fit using NNLO predictions. The fit is performed using the 
optimised renormalisation scheme for all energy points. 
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